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Dynamic  light  scattering  (or  photon  correlation  rnectroscopy.  PCS)  can  be  used  to  detenrune 
the  firequency  distnbution  of  the  light  scanered  from  polymer  solutions.  This  distnbuuon  cor. 
tains  information  about  the  dynamics  of  the  system.  Depending  on  the  wavelength  of  radiation 
and  the  size  of  the  polymer  molecule,  the  translations  and  rotations  of  the  enure  molecule,  or 
the  motions  of  the  monomers  can  be  studied.  This  information,  combined  with  static  light  scat¬ 
tering  (the  traditional  Zimm  Plot  analysis),  can  yield  a  molecular  weight  distnbuuon  tMWD), 
z-average  radius  of  gyration  (R^),  hydrodynamic  radius  (R|j).  translational  diffusion  cociTiciem 
(D,),  and  in  some  cases  a  rotational  diffusion  coefficient  (D^).  Dynamic  light  scattering  i  DLS  ; 
has  been  used  extensively  in  the  dilute  solution  regime  to  dcieimine  the  rv-pe  of  infomtation 
listed  above,  and  the  theories  associated  with  the  vanous  phenomena  are  weU  defined.  DLS  can 
also  be  applied  to  the  analyses  of  semi-dilute  soiuuons  and  to  study  the  dynamics  of  entangled 
systems  (e.g..  polymer  melts  or  diffusion  of  polymer  chains  in  a  matrix.),  although  the  interpre¬ 
tation  of  this  type  of  DLS  is  still  speculative.  This  report  will  focus  on  the  use  of  DLS  to  deter¬ 
mine  translational  diffusion  coefficients,  and  the  combination  of  static  and  dynamic  light  scat- 
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1.  Introduction 


1 .1  Historical  introduction 


The  initial  theory  of  light  scattering  by  gaseous  panicles  was  develojxjcl  m  ihe  hue  iX{K)\  b> 
Rayleigh,  and  extended  by  Mie,  Debye  and  others  in  the  early  1900’s  SmoluchovHski  <  Fklk  ’. 
and  Einstein  (1910)  developed  the  fluctuation  theory  of  light  scattering  which  accounted  for  the 
decreased  scattering  by  condensed  phases  due  to  destructive  interterencc  trom  scattering  ov  dil- 
ferent  molecules.^  The  application  of  light  scattering  to  detemiine  the  sizes  and  shapes  of  poh 
mers  c^me  to  maturity  m  the  194()’s  and  50’s.  uiitiaied  principally  by  the  work  ot  I>ehyc  and 
Zimm.^'  These  theories  all  dealt  with  the  intensity  of  scattered  light,  usually  measured  as  a 
function  of  a  number  of  scattering  angles,  i.e.,  'static"  light  scattering  Dynamic  light  scattering 
on  the  other  hand,  is  concerned  with  the  frequency  distribution  of  the  scattered  light  Tins  tre- 
quency  distribution  was  first  investigated  by  Brillouin  (1914  to  1920)  with  his  discovery  of  the 
inelastic  component  of  the  light  scattered  from  a  simple  liquid,  caused  by  scattering  from  ther¬ 
mal  sound  waves  in  the  fluid.  ^  This  "Raylcigh-Brillouin"  scattering  is  characterized  by  a  cen¬ 
tral  Rayleigh  peak  flanked  by  a  symmetric  Brillouin  doublet  '  (Figure  1 !.  where  the  frequency 
shift  of  the  Brillouin  peaks  (relative  to  the  central  Rayleigh  peak  i  is  proportional  to  the  velocity 
of  sound  in  the  fluid,  and  the  linewndth  of  the  Bnllouin  peak  is  related  to  the  attenuation  of 
sound  in  the  liquid.  Frequency  shifts  for  the  Brillouin  peak"^  are  on  the  order  of  0  ()1  w  avenum 
bers.  Raman  spectroscopy,  on  the  other  hand,  which  measures  the  frequency  shifts  of  the  scat¬ 
tered  light  caused  by  vibrational  rnottons  of  the  scattering  panicles,  exhibits  much  larger  shifts, 
typically  in  the  range  4000  -  400  cm'*.  The  linewidth  of  the  central  Rayleigh  peak,  which  is 
less  than  1  MHz  (10'^  is  much  smaller  than  either  the  Raman  or  the  Brillouin  shifts  Be¬ 

cause  the  frequencies  for  the  three  regimes  (Rayleigh  line-broadenmg.  Brillouin  shift.s  and 
linewidths,  and  Raman  shifts)  represent  such  a  large  energy  gap.  different  detection  schemes 
are  needed  for  each  technique.  Rayleigh  linewidths  are  generally  measured  with  a  digital 
correlator,  using  photon  countmg.  Both  the  frequency  shift  and  the  luiewidth  of  the  Brillouin 
peak  are  determined  with  a  Fabry-Perot  interferometer,  and  the  determuiation  of  Raman  fre¬ 
quency  shifts  is  usually  accomplished  via  a  grating  monochromator. 


Figure  1 .  Rayleigh-Briilouin  frequency  distribution  of  a  fluid 

Brillouin  spectroscopy  does  not  generally  have  much  application  in  polymer  solution  charac¬ 
terization,  since  the  breadth  of  the  Brillouin  peaks  is  very  narrow  in  polymer  solutions  Its  use 
in  polymer  analysis  is  generally  limited  to  light  scattering  of  bulk  polymers  (in  the  melt)  Tlie 
central  Rayleigh  peak  however  is  broadened  by  the  Brownian  motions  of  polymer  molecules  in 
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solution,  and  contains  dynamical  intormaiion  alx>ui  the  polviner  solution  Hie  obiccisM-  ot 
DLS  is  to  extract  information  on  the  translational  un  rotational)  motions  of  the  scattenne  me 
dium  from  the  linewidth  of  the  Rayleigh  scanered  hgnt  . 

1.2  Light  Scattering  Theory 

When  an  oscillating  electromagnetic  field  (such  as  in  a  beam  of  light)  approaches  a  non 
absorbing,  nonionizing  material,  an  oscillating  dipole  of  the  same  frequencs  as  the  nearb\  elec¬ 
tric  field  is  induced  in  the  material.  This  accelerating  dipole  then  radiates  energs  of  the  s,imc 
frequency  in  all  directions.  This  is  the  basis  for  the  phenomenon  of  light  scattering 

Figure  2  shows  a  schematic  view  of  a  typical  light  scattering  spectrometer^  (dyn,unit  or  static) 
Venically  polarized  laser  light  of  incident  frequency  passes  through  various  focusing  opnes 
and  impinges  on  a  small  scattering  center  in  a  sample  cell.  Tlie  scattered  light  at  an  angle  0 
(with  respect  to  the  exiting  beam  of  the  laser)  passes  through  two  pinholes  to  define  the  coher¬ 
ence  area,  and  is  collected  via  a  photomultiplier  tube  The  photon  pulses  are  amplified 

and  discriminated  and  then  fed  to  the  digital  correlator  for  analysis.  The  sample  cell  is  gener¬ 
ally  surrounded  by  a  vat  filled  with  a  liquid  whose  refractive  index  closely  matches  that  of  the 
sample  cell  (e.g.,  toluene)  to  limit  laser  flare  at  the  cell/air  interfaces. 

Scattering  of  electromagnetic  radiation  in  a  nonconducting,  nonmagnetic,  nonabsorbing  me¬ 
dium  occurs  as  a  result  of  local  fluctuations  in  the  dielectric  constant  tei  of  the  medium  'Hie 
dielectric  constant  is.  a  vector  quantity,  and  as  such,  consists  of  a  magnitude  and  a  directional 
component.  A  change  in  either  component  constitutes  a  change  in  the  dielectric  constant.  In  a 
pure  liquid,  the  magnitude  of  the  dielectric  constant  is  the  same  for  all  molecules,  but  the  direc¬ 
tions  of  the  dielectric  constants  keep  changing  due  to  the  Brownian  motions  of  the  molecules 
Thus,  scattering  in  a  pure  liquid  occurs  as  a  result  of  the  translational  motion  of  the  molecules 
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Figure  2.  Diagram  of  tight  scattering  apparatus 

1.2.1  Dynamic  Light  Scattering  (DLS) 

1. 2.1.1  Fluctuations  and  Autocorrelation  Functions 

Figure  3  shows  an  expanded  view  of  a  time-dependent  signal  A(t  ).  similar  to  what  the  scattered 
light  intensity  from  a  light  scattering  experiment  might  look  like:  The  apparent  randomness  of 
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the  signal  is  due  to  the  consnmi  motions  of  the  panicles  m  the  scattenne  region  Ov  er  a  long 
enough  time  period,  on  average  value  is  obtained,  but  m  the  microscopic  view,  there  appear  to 
be  random  fluctuations  in  intensity.  In  this  figure,  the  time  axis  has  been  broken  up  into  small 
sections,  :Sx  (small  compared  to  the  penod  of  the  fluctuations),  and  the  average  value  of  A, 
<A>,  is  equal  to  zero  (for  simplicity).  In  general,  the  values  of  the  signal  .4  at  two  differeni 
times,  Aj  and  A^^^  are  different,  but  if  the  tunc  increment  (_\t)  is  small  enough,  then  two  adjoin¬ 
ing  sections  of  the  signal,  Aj  and  Aj_^[.  will  have  values  close  to  each  other  Another  w  ay  of 
expressing  this  is  to  say  that  A.  and  A^^.^  are  "correlated"  over  a  small  time  penod.  but  that  the 
correlation  is  lost  as  T  approaches  and  exceeds  the  time  scale  of  the  fluctuations.  The  length  of 
time  over  which  the  signal  is  correlated  is  mherent  in  the  data.  The  autocorrelation  funciion  oi 
a  signal.  A(t).  is  expressed  mathematically  by 

C(t)  =  (A(0)A(r))  =  |.j^A(t)/A(r  +  r)di ,  <!> 

where  dt  is  approximately  equivalent  to  the  discrete  interval  At  described  previously  Tins 
autocorrelation  definition  is  the  description  of  one  point  (X)  on  the  autocorrelation  function  -  it 
is  the  sum  of  a  series  of  multiplications.  The  entire  autocorrelation  function  is  constructed  by 
computing  <A(0)A(t)>,  for  z  ranging  from  0  to  a  point  at  which  the  sign,il  is  no  longer 
correlated. 


Figure  3.  Time-dependent  property  "A(t)" 

Some  of  the  points  in  Figure  3  are  positive,  and  some  are  negative;  thus  the  various  products 
computed  in  Equation  I  will  be  positive  as  well  as  negative  (e  g..  .Aj  x  A,^  is  negative,  but  x 
Aj,  and  A^  x  A^  are  positive).  The  first  point  of  the  autocorrelation  funciion  is  <,4(()).4(0)>. 
where  Z  in  Equation  I  is  replaced  by  0.  This  point  is  actually  the  sum  of  the  square  of  every 
point  along  the  curve.  Since  the  square  of  a  number  is  always  positive,  the  first  point  of  the 
correlation  function,  <A(0)A(0)>,  wUl  be  the  largest  value  in  the  correlation  function,  because 
it  will  not  contain  any  negative  components. 

The  point  <A(0)A(1)>  can  be  calculated  by  shifting  the  A(t)  curve  one  unit  ( 1  At)  to  the  nght. 
and  computing  the  sum  of  the  products  of  where  the  two  curves  overlap.  The  curv'es  are  nearly 
superimposed,  so  the  value  A(0)A(1)  will  be  only  slightly  less  than  A(0)A(nj.  since  there  are 
only  a  few  negative  products.  Eventually,  as  x  gets  large,  there  is  very  little  overlap  between 
the  two  curves,  and  an  averaging  process  takes  place,  luid  the  correlation  function  would  ap¬ 
proach  zero  (or  <A>'',  if  <A>  0). 


3 


The  actual  correlation  function  that  would  be  computed  from  this  hypothetisa]  Mcnal  An*,  u  iH 
look  essentially  like  an  exponential  decay  curv'e  (Figure  4  shows  an  autocorrelation  function 
from  a  light  scattering  experiment.)  If  the  correlation  function  can  tse  described  as  a  sirndc 
exponential  decay,  (i.e.,oniy  one  physical  proce.ss  is  causing  the  fluctuations!.  Equation  I  could 
be  simplified  to 


(A(0l4(r))  =  (A)  exp  ~  <  -» 

where  is  the  characteristic  decay  tune  associated  with  whatever  phenomenon  is  causing  tfve 
fluctuations  in  the  signal  Ad).  If  the  panicles  diffusmg  in  the  medium  were  all  of  the  same 
size,  they  would  have  the  same  diffusion  coefficiem.  and  a  single  decay  time,  t i  related  to  the 
diffusion  coefficient  of  the  molecules)  could  be  extracted  from  the  curve  How'ever,  since  mac- 
romoiecules  always  have  some  degree  of  poly dispers tty.  there  will  be  more  than  one  diffusion 
coefficient  contributing  to  the  autocorrelation  function,  and  the  curve  actually  consists  of  a 
superposition  of  many  exponential  decay  curves. 

230-1 
C(t)  221  - 

Arbitrary 
Units  ' 

203- 

1954 

0  28  56  84  112  140 
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.Autocorrelation  function  measured  for  sample  BR3M  (cone.  =  0.5035 
mg/mL)  at  scattering  angle  of  30°  and  sample  time  ( Ax)  =  48  psec. 

Figure  4.  Autocorrelation  function  of  scattered  light  intensity 

Since  Equation  2  is  a  function  of  time,  the  Fourier  transform  of  this  equation  should  produce  m 
equation  in  the  frequency  domain,  termed  the  spectral  density,  or  power  spectrum.  This  is  the 
spectrum  of  the  light  scattered  by  the  molecules  in  solution,  as  shown  in  Figure  I,  and  the 
linewidth  (F  =  2ji/x^),  can  be  extracted  from  the  experimentally  determined  x_,.  .As  in  the  corre¬ 
lation  function,  this  Lawrencian  curve  is  actually  a  superposition  of  many  Lawrencian  curs  es 
with  different  linewidths,,  all  centered  at  the  incident  frequency.  The  extraction  of  multiple 
linewidth  from  an  experimental  correlation  function  will  be  discussed  in  Section  1 .3. 1 

1 .2.1 .2  Diffusion  Coefficients  from  Autocorrelation  Functions 

The  scattered  light  intensity.  from  a  solution  of  macromolecuies  impinges  on  a 

photomultiplier  tube  which  in  turn  outputs  a  current,  i(t).  proportional  to  I^.(t).  Electromagnetic- 
theory  tells  us  that  the  light  intensity  is  equivalent  to  the  square  of  the  electric  field;  thus  it  fol¬ 
lows  that  I^(t)  =  i(t).  The  output  of  the  PMT  is  fed  into  a  correlator  which  computes 
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the  intensify  autocorrelation  function. 


(i{0)i(r))  =  B  ( i£s(0i  i  ■ 

where  B  is  a  proponionalirv  constant.  Equation  3  is  a  fourth  degree  correlation  hinctiun.  u  Inch 
is  difficult  to  analyze.  In  the  case  of  dilute  ptslymer  solutions,  however.  follow  s  a  Gaussian 
distribution,  and  the  correlation  function  can  be  simplified  to  the  sum  of  two  correlation  func¬ 
tions,^  the  first  of  which  is  a  constant  term.  <E^(0)E^(0)>'  (the  baselme  of  the  correlation  func 

tion),  and  the  second  term  which  is  a  simpler  second  destree  correlation  function. 
*  *  ■> 

<E^  {0)E^  (t)>".  the  electric  field  autocorrelation  function. 

Incident  light  which  impinges  on  a  molecule  w  ith  a  polarizability  tensor  a  induces  a  tluctuaiuic 
dipole  moment.  The  scattered  electric  field  at  the  detector  ts  proportional  to  te^^ 
where  ct  Jt)  is  the  component  of  the  polarizability  along  the  initial  and  final  polanzation  direc¬ 
tions,  ana  the  factor  varies  when  the  molecule  translates,  where  rtt)  specifies  the  mole¬ 

cule’s  position  at  time  t.  The  wave  vector,  q.  is  equal  to  (47in//jsint9/2).  where  n  is  the  refrac¬ 
tive  index  of  the  medium.  X  is  the  wavelength  of  the  incident  radiation,  and  0  is  the  scattering 
angle.  Thus  the  electric  field  autocorrelation  function.  <£,,  (O)E,^  (t)>.  is  related  to  the  molecu¬ 
lar  polarizability. 

When  the  product  of  the  wave  vector  and  the  radius  of  the  scattering  panicle  is  less  th;ut  unity 
(qR  «  1).  the  tune  correlation  functions  are  sensitive  to  fluctuations  occurring  on  a  nme-scaJe 
associated  with  center-of-mass  diffusion.^  When  q  is  increased  (by  decreasing  the  w  aselength 
or  increasing  the  angle)  and  qR  exceeds  1,  then  internal  motions  'rotations)  of  the  polymer 
molecules  become  imponant.  This  is  an  intermediate  q  range,  and  the  interpretation  of  DLS 
data  as  due  strictly  to  center-of-mass  diffusion  can  be  erroneous.  Smce  the  wavelength  of  light 
cannot  be  readily  changed.  DLS  measurements  of  large  particles  must  be  made  at  low  iingles  to 
exclude  any  rotational  contributions  to  the  correlation  function.  Finally,  when  q  gets  very  large 
such  that  qa  is  on  the  order  of  1,  where  a  is  the  size  of  polymer  repeat  unit,  then  the  time  corre¬ 
lation  functions  are  probing  motions  of  the  monomers.  This  requires  wavelengths  on  the  order 
of  a  few  Angstroms,  the  frequency  of  neutron  radiation,  and  thus  not  accessible  through  light 
scattering  with  visible  radiation. 

If  one  can  assume  that  the  experimentally  measured  time-correlation  function  is  caused  by 
center-of-mass  diffusion,  then  the  diffusion  coefficient  of  the  polymer  molecule  can  be  simply 
related  to  the  linewidih  via  T  =  The  diffusion  constant  thus  deienmined  can  be  related  to 

the  size  of  the  diffusing  panicle  through  the  Stokes-Einstein  equation. 

l^gT 

i4) 

h 

where  kg  is  the  Boltzmann  constant.  T  is  the  temperature,  q  is  the  viscosity  of  the  medium  the 
panicle  is  diffusing  in,  and  Rj^  is  the  hydrodynamic  radius  of  the  panicle.  For  sphencal  pani¬ 
cles  {e.g.  latex  spheres)  DLS  can  be  used  to  measure  actual  panicle  sizes,  however,  for  dilute 
polymer  solutions,  this  is  really  an  "equivalent"  hydrodynamic  ladius,  with  no  physical 
significance. 


1.2.1 .3  Molecular  Weight  Distributions  from  Diffusion  Coefficients 

Since  the  time  correlation  function  measured  from  a  DLS  exfseriment  contains  miormauon  on 
many  difterent  linewidths  (due  to  the  polydispersity  of  the  polymer  sample),  abstracuon  oi  a 
distribution  of  linewidths  representative  of  the  different  polymer  molecules  m  the  sample 
should  be  possible.  This  distribution  of  linewidths,  GlP),  can  be  transfonned  to  a  distribution 
of  diffusion  coefficients,  G(D),  which  then  can  be  transfonned”  to  a  distribution  of  niolecuLu 
weights.  f(M). 

The  area  under  the  distribution  curve  GtT)  is  proportional  to  the  tune-averaeed.  total  iniensits 
of  light  scattered  at  an  angle  0  (deieimined  from  a  static  light  scattering  e.xpenment  >. 

OO 

(fio))  =  j  G(r>dr 

This  time-averaged  intensity  can  also  be  expressed  in  terms  of  the  molecular  w. eight  compo¬ 
nents  of  the  polymer  species  by  the  following  equation^  from  classical  light  scattering  analysts. 

where  K‘  contains  a  collection  of  optical  constants.  C|  is  the  concentration  i  mass/vciiume  >  of  a 
polymer  molecule  of  molecular  weight  M^,  and  P(0.M|  ,i  is  the  panicle  scattering  factor  vs  htch  ts 
different  for  different  shaped  molecules  (i.e..  rigid  rods,  random  coils,  etc.). 

If  C-  can  be  equated  to  f(M')AM'/V,  where  f(M-)  represents  a  molecular  weight  distribution 
curve,  AMj  is  a  small  increment  along  the  curve  and  V  is  a  unit  volume,  then  the  distribution  of 
linewidths  can  be  related  to  the  distribution  of  molecular  weights  via 


This  transformation  from  gamma-space  to  molecular  weight-space  rakes  place  by  equating 
small  areas  under  the  Gif)  curve.  AT,  with  area  under  the  f(M)  curse,  AM  The  diffusion  coef¬ 
ficient  has  a  concentration  dependence  which  can  be  expressed  as 

JDiC)  ~  Z?(j(  1  +  kpC)  (k) 


where  kp  is  the  second  virial  coefficient  for  diffusion.  can  be  estimated  by  extrapi^laimg 
measured  diffusion  coefficients  at  various  concentrations  to  infinite  dilution.  The  diffusion  co¬ 
efficients  should  also  be  measured  at  several  angles,  particularly  for  high  MW  polymers,  and 
extrapolated  to  zero-angle  to  exclude  any  rotational  contributions  at  higher  luigles.  T!ie  infinite 
dilution  diffusion  coefficient  can  then  be  related  to  a  molecular  weight  via 

Do  =  .  (9) 


where  the  pre-exponential  factor,  kj,  depends  on  the  polymer/solvent  system,  and  must  be  de¬ 
termined  experimentally.  The  exponent  b  is  equal  to  0.5  for  random  coils  at  the  theta  contiinon 
and  has  been  verified  experimentally.  For  polymers  in  good  solvents  h  is  expected  to  he  closer 
to  0.6,  although  experimentally  reported  values  vary  from  0.55  to  0.68.^'^ 
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If  kp,  kj  and  b  ;ire  known  (or  estimated)  then  molerular  weights  can  be  calculated  ttom  expen- 
mentally  determined  ILnew  idths  via  the  following  equation 

M  =  a-,, 

where  the  values  of  M  can  be  substituted  into  Equation  7  for  the  N1  values  and  the  AM  mere 
ments  can  be  calculated  by  determining  the  values  of  M  at  the  endpoints  of  the  AF  increments 

Each  M|  in  Equation  7  must  have  a  corresponding  which  is  used  to  detemiine  Pf8,  and 
so  a  radius  of  gyration  is  computed  for  each  moleculim  weight  similarly  to  Equation  '■K 


Like  b.  the  exponent  b’  is  expected  to  be  around  0.6  for  random. i  coils  in  a  good  solvent,  but  nut 
necessarily  equivalent  to  the  dynamic  expionent  b.  The  pre-exponential  kj‘  depends  on  the  par¬ 
ticular  chemical  system.  For  random  coil  molecules,  the  panicle  scattering  factor  usually  takes 
the  form 

PUC.M)  =  X  le-^-  1  +  .rj  r/:, 

X- 

where  X  =  q^R^j",  and  R^  is  the  radius  of  gyration  computed  from  Equation  1 !  B\  combining 
the  results  of  static  and  dynamic  light  scattering  measurements  with  estimates  of  kj.  k’y,  b.  and 
b',  an  iterative  process  can  be  used  to  determine  an  effective  molecular  weight  distribution. 

As  can  be  seen  from  the  preceding  discussion,  much  information  about  the  polymer-solvent 
system  must  be  known  in  order  to  determine  a  molecular  weight  distnbution  in  this  manner 
Additionally,  most  of  the  equations  apply  strictly  to  monodisperse  species,  and  any  breadth  to 
the  molecular  weight  distribution  adds  complications.  For  these  reasons,  this  method  of  deter¬ 
mining  a  molecular  weight  distribution  is  not  of  general  utility  but  in  certain  instances,  may  be 
appropriate. 

Size  exclusion  chromatography  (SEC)  is  an  experimentally  simpler  method  of  detemiining  mo¬ 
lecular  weight  distributions,  but  has  several  limitations.  SEC  requires  the  use  of  polymer  stan¬ 
dards  to  determine  actual  molecular  weight  values,  and  well-charactenzed  standards  are  onl\ 
available  for  very  common  polymers.  Light  scattering,  however,  is  an  absolute  technique  and 
requues  no  standards.  Additionally,  certain  polymers  may  only  be  soluble  in  solvents  which 
•  are  incompatible  with  the  chromatography  coluirins.  or  at  temperatures  not  accessible  by  ordi¬ 

nary  SEC  instrumentation.  For  example,  aromatic  polyamides  such  as  Kevlar,  and  polycuy  ! 
ether  ketones)  are  generally  .soluble  only  in  strong  acids.  Corrosive  solvents  are  not  a  problem 
for  light  scattering  which  is  usually  measured  in  a  glass  or  quartz  cell. 

1.2.2  Static  Light  Scattering  (SLS) 

q 

Classical  electromagnetic  theory  and  solution  thermodynamics  show  that  the  intensity  of  light 
scattered  (from  a  polarized  light  source  of  wavelength  X)  by  small  (compared  to  A),  isotropic 
scatterers  at  an  angle  0  can  be  represented  by 


7 


The  Rayleigh  ratio.  R^,  contains  the  fX'ntnent  measured  quantities,  .uui  is  equal  iti  i^^r"  I  .  '.uierc 
Irt  i-S  the  measured  scattered  light  tritensity.  the  meuieni  iniensii>,  and  v  is  ihe  c^t^>.cn!ra!ios| 
of  the  solution  i mass/ volume  i  Other  parameters  are  grouped  unto  the  opuc.il  constant  K.  sic 
fined  as  4?t''n"(da'dc)*/(LA'*i.  where  n  is  ihe  refractive  inde.s  of  the  medium,  dii/dc  sthe  tefiat 
live  tnde.x  increment)  is  the  change  in  refractive  mde.\  of  the  solutuui  svnh  concentraiion.  .uui  L 
IS  Avogadro’s  number,  dn  practice,  the  scattering  distance  r  and  the  incideni  intensity  I  ^  are 
not  measured  Instead,  a  solvent  for  which  this  Rayleigh  ratio  has  been  picviousls  Jetennined 
{toluene  or  benzene,  rypically!  is  used  as  a  calibration  liquid  Its  scatierc  ‘  miensits  is  mea¬ 
sured.  and  the  known  value  of  us  Rayleigh  ratio  is  used  to  calculate  ihc  v,ilucs  o'  t  and  !  _  hu 
the  current  system,  > 

Equation  13  shows  the  concentration  dependence  ot  the  scattered  light  mtcnsits  ot  .i  [xihn.cr 
.iolution.  For  large  (r  >  1/20  A)  molecules,  there  is  also  an  angular  dependence  of  the  scattered 
mtensiry  due  to  the  size  of  the  molecule  Small  scattcrers  have  a  symmetric  scattering  pattern 
about  0  =  90^.  because  they  look  like  a  pomi  to  the  incident  radiation  L.uge  particles,  ho\s- 
ever.  will  scatter  radiation  from  different  ponions  of  their  molecule,  and  thus  there  will  lx  de¬ 
structive  interference  at  the  detector  due  to  the  phase  differences  of  light  arriving  at  the  detector 
which  has  scattered  from  two  different  points  on  the  same  molecule  The  scattered  light  inten¬ 
sity  will  be  highest  at  large  angles,  because  the  difference  m  path  length  is  larger,  .uiii  this 
assyrrumerrs'  about  dO'^^is  used  to  extract  size  infomiation 


The  angular  dependence  of  the  scattered  light  intensity  takes  the  tollowmg  fomi  -  the  ellipsis 
indicates  that  higher  order  tenns  m  R^'sin'(0/2)  arc  present 


Kc  _ 


'vi* 


For  high  molecular  weight  polymers  with  large  radii,  the  curvature  in  Equation  U  can  be  sig¬ 
nificant.  and  thus  the  initial  slope  must  be  used  to  obtain  accurate  values  of  ,\f^  and  liqua¬ 
tions  13  and  14  can  be  combined  and  the  molecular  weicht.  second  viriai  coefficieni'  .uid  the 


radius  of  gyration  can  be  obtained  from  a  single  plot,  to  be  described  in  section  *  3  2  The  mo¬ 
lecular  weight  determined  is  absolute  (i.e  .  independent  of  the  solvent  used),  and  is  the  eight- 
average  value,  as  oppo.sed  to  the  number-average  .M^  determined  from  colligative  propenv 
measurements.  The  second  virial  coefficient  is  a  measure  of  the  thermodynamic  interaction  ot 


the  solute  and  solvent  lihc  third  virial  coefficient  is  usually  small  and  not  impon.int  unless  the 
solution  concentraiion  is  high).  The  radius  of  gyration  determined  from  light  scattering  mca 
surements.  defined  as  the  root  mean  square  distance  of  the  segments  of  the  molecule  from  iis 
center  of  mass,  is  a  z-average  value  The  size  of  the  molecule  and  the  second  viruil  coefficient 
are  not  absolute  values,  but  depend  on  the  solvent  used.  .A-,  is  positive  for  a  thennodynarnicaily 
good  solvent  (polymer-solvent  affinity  is  greater  than  polymer-polymer  affinity  >.  ai  i  s  anes  m 
magnitude  for  different  solvent  qualities,  is  negative  for  a  poor  solvent,  and  zero  (equal  affini¬ 
ties)  at  the  theta  condition.  Accorduiglv.  the  radius  of  gyration  is  largest  in  a  good  sohent. 
smallest  m  a  poor  solvent,  and  an  intermediate  value  in  a  theta  solvent. 
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1 .3  Data  Analysts 

1.3.1  Dynamic  Light  Scattering 

In  dynamic  light  scanermg  experunents.  there  are  two  autcKorrelation  funcMofis  ot  unerest,  the 
intensity  autocorrelation  function  g'*’(t).  which  is  the  function  actually  computed  by  the 
correlator  (Equation  3),  and  the  electric  field  autocorrelation  function,  g  '  (ti.  which  .s  related  to 
the  molecular  parameters,  desenbed  in  section  1.2  1.2.  These  two  correlation  functions  axe  re¬ 
lated  through  the  Siegen  relation, 

=  1  ^  {S 


where  the  baseline  of  the  experimental  correlation  function  (Ai  has  been  subtracted  The  con¬ 
stant  p,  0  <  p  <  1,  is  usually  determined  as  a  parameter  during  the  fitting  of  the  data  Tlic  eicc- 
tric  field  correlation  function,  g‘''(t),  is  assumed  to  be  in  the  fomi  of  a  sum  of  singk- 
exponentials, 

g"'(r)=|^  GiDe-^^dT 


where  the  T’s  represent  the  different  Imewidths  due  to  the  different  molecular  weight  sfwcies 
and  G(n  represents  theu  distribution  function.  As  mentioned  in  sections  i  2  I  2.  in  the  ab¬ 
sence  of  internal  motions,  the  diffusion  coefficient  is  related  to  the  linew  idth  via 


r  -  Dj-q-.  i  }7) 

The  goal  of  any  light  scattering  data  analysis  routine,  therefore,  i-s  to  inven  Equation  16  to  ob¬ 
tain  There  are  many  mathematical  methods  available  for  this  transformation  -  a  few 

methods  which  are  available  on  the  Brookhaven  instrument  will  be  reviewed  brieflv  ,\J1  of 
the  methods  described  here  are  least-squares  fits,  m  which  the  objective  function  to  I-k:  mini¬ 
mized  is. 


53  -  y*(MT)]- 


/=! 


,!S! 


where  y^jdAt)  is  the  measured  correlation  function  (with  the  ba,seline  A  subtracted)  and 
y^flAT)  is  the  model  correlation  function  proposed  by  the  panicular  method  of  data  analysis 
The  index  I  mns  over  all  the  data  points  in  the  coneiation  function  tin  the  Brookhaven 
correlator,  there  are  128  data  channels).  Equation  18  is  minimized  with  respect  to  a  number  of 
parameters  aj  (j  =  1.  2 . M),  which  vary  for  the  different  methods  of  data  analysis 

1 .3.1 .1  Method  of  Cumulants 

The  method  of  cumulants' ’  is  the  simplest  way  to  analyze  DLS  data,  but  also  yiekls  the  least 
amount  of  information.  Only  the  average  ILnewidth  and  its  variance  can  be  obtained  tthe  first 
and  second  cumulants)  with  any  degree  of  accuracy.  (Sometimes  the  3rd  cumulant  is  also  re¬ 
ported.)  The  cumulants  are  defined  in  terms  of  moments  about  the  mean  linewidth.  where  the 
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moment  is  defined  as 


u 


m 


chFnr-  r  I  dT 


and  the  resulting  computed  correlation  function  is 

y  *  (r)  =  (A^)^  'exp (^- Ft  +  4^2^' 


f/M) 


(2U) 


where  the  ellipsis  indicates  that  higher  cumulants  may  be  included  The  parametcrN  used  to 
minimue  Equation  1 8  are  the  cumulants  (F,  p-,,  etc.)  and  the  factor  Ap 


1.3.1 .2  Exponential  Sampling 

Another  common  data  analysis  routine  used  in  DLS  is  the  Exponential  Sampling  Technique, 
also  called  the  LaPlace  Transform  Inversion.*'  which  is  based  on  the  eigenfunctions 
eigenvalues  of  the  LaPlace  transform  of  G(  F), 

\p^(F)exp(-Fx)dr  (2!  I 

where  analytical  expressions  for  and  are  given  m  Reference  1 2 

The  distribution  function  GiF)  is  then  expanded  into  its  complete  set  of  eigenfunctions  via 


G{F) 


y.<^{F)  dio 


By  substitution  of  Equation  22  into  Equation  16  and  using  the  relationship  in  Equation  21.  the 
correlation  function  can  be  written  as 


a^^xp^ix)  dw 


(2Si 


As  (o  gets  very  large  (to  >  the  eigenvalues  approach  zero,  and  their  contribution  to  the 

correlation  function  cannot  be  distinguished  from  experimental  error,  thus  the  distribution  of 
linewidths  is  band-limited  by  changing  the  limits  in  Equation  22  from  ±=0  to  After  a 

series  of  mathematical  manipulations,  GiF)  may  be  represented  as  a  series  of  exponentially 
spaced  samples  in  F-space,  with  amplitudes  a^^.  The  amplitudes  are  used  to  minimize  Equation 
20.  with  the  proposed  correlation  function  given  by 

M 


y  *  (Mr) 


a„  ex{A-FmlAx) 


i24) 


The  best  fits  to  the  data  are  usually  obtained  with  greater  than  about  3. 

1.3.1 .3  Non-negatively  Constrained  Least  Squares 

This  data  analysis  routine  is  an  adaptation  of  the  previous  method,  developed  for  multimodal 
linewidth  distributions.*^  Equation  24  is  used,  and  each  a^^  is  constrained  to  be  non-negative.  In 
this  method,  the  smallest  and  the  largest  linewidths  (r_.,„  and  F  )  are  predefined  to  be 
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functions  of  the  experimental  sample  tune  and  the  average  luiew  idih  dcicnnuicd  trurn  the  dojsc 
of  the  loganthm  of  the  measured  data,  in  the  Broukhaven  version  of  this  data  rouunc. 

the  user  is  given  the  option  of  overriding  the  default  and  values 

1.3.1 .4  Double  Exponential 

This  method  forces  the  correlation  function  to  be  composed  of  two  discreet  exponential  func¬ 
tions,*'*  as  shown  in  the  following  equation. 

=  /1|  expf-r^r)  +.4-,  expi-An. 

This  is  a  dangerous  method  for  analyzing  correlation  functions  unless  the  user  is  cenain  of  the 
bimodal  character  of  the  data,  because  any  noise  in  the  data  could  be  mterpreted  as  anuilicr  spe¬ 
cies.  This  type  of  data  analysis  is  primarily  used  m  DLS  of  entangled  systems,  where  a  fast 
cooperative  diffusion  mode  of  a  network,  and  a  slower  self-diffusion  mode  of  a  smaller  unit  is 
present.*^ 

1.3.2  Static  Light  Scattering 

For  dilute  solutions,  the  measured  scattered  light  intensity  from  a  given  illuminated  area  (the 
scanering  volume)  is  the  superposition  of  the  scanered  light  intensities  from  all  elements  within 
that  region.  Thus,  the  scattered  light  measured  from  a  solution  of  polymer  molecules  consists 
of  the  light  scattered  by  the  small  solvent  molecules  plus  the  scanering  from  the  large  poly¬ 
meric  species.  If  the  scattering  of  the  pure  solvent  is  measured  separately,  it  can  be  subtracted 
from  the  solution  scattering  and  the  scanering  due  to  the  polymer  molecules  alone  can  be  iso¬ 
lated.  The  resulting  excess  scattered  intensities  are  convened  to  Rayleigh  rafio.s.  as  described  in 
section  1.2.2.  The  data  is  usuaUy  analyzed  via  a  Zimm  Plot  in  which  the  (KC/R^)  values  are 
plotted  on  the  y-axis  versus  (sin"(0/2)  +  kc)  on  the  x-axis,  where  k  is  an  arbitrary  constant  used 
to  spread  out  the  points  on  the  plot  The  set  of  points  corresponding  to  the  same  concentration 
(different  scanering  angles)  should  form  a  straight  line  parallel  to  the  other  concentrations.  An¬ 
other  set  of  parallel  lines  (usually  with  smaller  slopes)  can  be  drawn  by  connectuig  the  points 
corresponding  to  a  given  scattering  angle  at  different  solution  concentrations.  These  two  sets  of 
parallel  lines  form  a  grid-like  structure  (see  Figure  5). 

The  constant  concentration  lines  can  be  extrapolated  to  a  "zero-angle"  pomt  (with  an  x- 
coordinale  equcil  to  kc).  Likewise,  the  constant  angle  lines  can  be  extrapolated  to  a  "zero- 
concentration"  point  (with  an  x-coordinate  equal  to  (sin"(0/2)).  Each  of  the  extrapolated  "zero- 
angle"  points  can  be  extrapolated  to  zero-concentration  (x  value  equal  to  0);  similarly  the  ex¬ 
trapolated  "zero-concentration"  points  can  be  extrapolated  to  zero  angle  (x  value  equal  to  0). 
and  these  two  doubly  extrapolated  lines  should  intersect  at  the  same  pomt  on  the  y-axis.  the  re¬ 
ciprocal  of  the  molecular  weight.  The  second  virial  coefficient  is  obtained  from  the  initial  slope 
of  the  zero-angle  line  (see  Equation  13)  and  the  radius  of  gyration  is  abstracted  from  the  initial 
slope  of  the  zero-concentration  line  (Equation  14).  The  initial  slopes  of  both  of  these  lines  must 
be  used,  as  higher  order  terms  become  significant  as  concentration  and  angle  increase.  The  cur¬ 
vature  of  the  zero-angle  line  should  be  negligible,  unless  the  solutions  are  too  concentrated,  hut 
the  curvature  of  the  zero-concentration  line  can  be  significant,  particularly  for  large  polymeric 
species.  A  variation  of  the  Zimm  Plot  used  for  analyzing  high  molecular  weight  samples,  was 
proposed  by  Berry, *^  is  to  plot  (KC/Rg)*^“  values  on  the  ordinate  vs.  the  same  abscissa  values. 
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The  intercept  then  yields  the  reciprocal  of  (M^,r'“.  and  ih  *  equations  for  detemumni:  A-,  and 
Rj,  are  slightly  different. 

The  larger  the  difference  in  scattered  intensity  between  the  solvent  and  solution,  the  greater  the 
precision  of  the  parameters  determined  from  the  light  scattering  experiment  The  increase  in 
scattering  intensity  of  the  polymer  solution  over  that  of  the  pure  solvent  (at  constant  measuring 
angle)  depends  on  four  factors;  (1)  the  polymer  molecular  weight,  (2)  the  concentration  of  the 
polymer  solution.  (3)  the  difference  m  the  refractive  index  between  the  solvent  and  the  pKilymcr 
solution,  and  (4)  the  wavelength  of  the  laser.  The  piolymer  molecular  weight  is  obviously  fixed, 
but  the  other  parameters  can  be  adjusted  to  increase  the  scattering  intensity  for  a  poor  scatterer 
Increasing  the  solution  concentration  (bemg  careful  to  remain  in  the  dilute  solution  regime  i. 
changing  the  solvent  to  one  with  a  larger  dn/dc  value  for  the  panicular  prolyrner.  and  decreasing 
the  wavelength  of  the  laser,  all  increase  the  scattermg  intensity. 


2.  Experimental 
2.1  Sample  Preparation 


Table  1.  Polymer  samples  specifications 


MwxlO® 

. "*1 

Mw/Mn 

Composition 

N170K 

0.168 

1.04 

100%  PL  #20 137-2  I 

N770K 

0.765 

1.04 

100%  PL  #20140-9  1 

0.990 

1.04 

100%  PL  #20141-5  ; 

1.53 

1.06 

100%  PL  #20142-5  1 

N2M 

2.14 

1.06 

100%  PL  #20143-2 

2.28 

1.05 

100%  PL  #20 143 -7  1 

N3M 

2.91 

1.04 

100%  PL  #20145-9 

4.06 

1.06 

100%  PL  #20146-6  I 

9.35 

1.20 

100%  PL  #20148-3  1 

BR2M 

1.85** 

I  2* 

47%  PL  #20142-5,  53%  PL  #20143-2 

BR3M 

3.17^* 

1.2* 

56%  PL  #20143-7  44%  PL  #20146-6  I 

BI850K 

0.839** 

2.0* 

66%  PL  #20137-2,  34%  PL  #20143-2  j 

BI4M 

4.08** 

2.0* 

63%  PL  #20 14 1-5,  37%  PL#20 148-3  j 

♦Estimated  values  of  and  | 

i 

.  .  . .  ,  _  _  . . . . .  J 

Polystyrenes  (PS)  with  three  types  of  molecular  weight  distributions  (narrow,  broad,  and 
bimodal)  were  investigated  to  demonstrate  the  capabilities  of  dynamic  and  static  light 
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scattering.  Several  commercial  PS  standards  from  Polymer  Laboratories  vs  etc  used  auiivuiu- 
ally,  and  in  various  combmattons  and  then  propenies  are  listed  ui  l  ablc  i  bout  viaiuianls 
(N170K.  N770K,  N2M.  and  N3M)  were  used  xs  is  for  the  narrow  MW  dtsirubnon  samples 
Two  broad  distributions  (BR2M  and  BR3M)  and  two  bimodal  distnburtons  (B1850K  .md 
BI4M)  were  simulated  by  mixing  PS  standards  with  similar  and  dissumlar  molecular  weights, 
respectively. 

Stock  solutions  of  concentration  approximately  1.0  mg/mL  of  polystyrenes  were  prepared  w  ith 
toluene  (EM  Science,  spectroscopic  grade),  and  allowed  to  dissolve  at  room  temperature,  vvith 
slight  agitation,  for  at  least  24  hours.  These  stock  solutions  were  diluted  to  three  or  four  con¬ 
centrations  of  0.1  -  0.5  mg/mL  for  static  light  scattering  experiments,  except  for  the  lowest  mo¬ 
lecular  weight  samples  in  which  the  highest  concentrations  were  about  10  mg/mL  Concentr.]- 
tions  used  in  dynamic  light  scanenng  were  slightly  higher  -  about  0.5  to  1  {)  ing/mL 

2.2  Instrumentation 

All  light  scattering  measurements  were  made  on  a  Brookhaven  Instruments  BI2(K)SM  goniome¬ 
ter  and  BI2030AT  128  channel  correlator  using  a  Melles  Griot  5  mW  helium-neon  laser  (6.^3 
nm)  as  the  light  source.  The  refractive  index  increment  (dn/dc)  of  PS  in  toluene  was  not  mea¬ 
sured.  but  a  literature  value*  ‘  of  0.107  mL/g  was  used.  Toluene  was  also  used  as  the  calibration 

4  1  O  ^  t 

liquid,  with  a  Rayleigh  ratio  of  1.41  x  10''  cm’  .  Solvent  and  polymer  .solutions  were  filtered 
through  0.2p  and  0.45p  PTFE  filters,  respectively.  The  temperature  of  the  light  scattering  cell 
was  maintained  at  25.0  ±0.1  C  for  all  measurements. 


3.  Results  and  Discussion 
3.1  Static  Light  Scattering 

Static  light  scattering  measurements  were  made  at  angles  30°,  37.5^’,  45^.  60”.  75'’,  90'’.  1 05'’. 
120®,  135®,  142.5®,  and  150®  on  the  eight  polymer  solutions  described  ui  the  previous  section. 
The  static  LS  parameters  are  listed  in  Table  2  along  with  the  data  analysis  method  used. 

Figure  5  shows  typical  Zimm  and  Berry  plots  obtained  for  different  polystyrene  samples.  The 
lower  molecular  weight  samples  showed  a  linear  concentration  dependence  and  thus  were  best 
analyzed  through  a  classical  Zimm  plot  approach,  however,  a  curvature  to  the  fixed  concentra¬ 
tion  data  was  observed  for  >  1.0  x  10®,  and  thus  a  Berry  plot  gave  a  better  fit  to  the  data. 
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Table  2.  Static  light  scattering  parameters 


g/  mol 

R  ,  nm 

9 

10“"  X  A., 

2 

ml  mol^g 

Type  of  Plot’ 

0.178 

Zunm 

0.168 

Debye 

N770K 

0.73 1 

36.6 

4.(K) 

Zimm 

BI850K 

0-729 

63.4 

0,537 

Zmun 

BR2M 

1.94 

2.52 

Berry 

N2M 

2.21 

74.7 

2.44 

Berry 

BR3M 

3.01 

91.6 

. 

2.14 

i  Berry 

N3M 

3.07 

89.8 

2.24 

;  Berry 

BI4M 

3.49 

129.8 

1.83 

1  Berry 

*Zimm  Plot:  linear  fit  to  concentration  and  angular  data; 

Debye  Plot:  linear  fit  to  90°  data  only; 

Berry  Plot:  second  degree  fit  to  concentration,  linear  fit  to  angular  data. 
**Less  than  the  theoretical  lower  lunit  of  detection  of  =  a.^/20  =  32  nm 


Zimm  Plot  (left)  for  N770K  (cone.  =  0.2505.  0.5010,  and  0.7515  mg/mL);  and  Berry' 
Plot  (right)  for  N3M  (cone.  =  0.1007,  0.2014,  0.3021,  and  0.5035  mg/mL). 


Figure  5.  Typical  Zimm  and  Berry  Plots 

.  .  The  Rg  and  A2  values  for  the  unimodal  PS  samples  were  linear  with  molecular  weight  in  log- 
log  plots  as  can  be  seen  from  Figure  6.  The  relationships  derived  from  these  log-log  plots  are 
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expressed  mathematically  by  the  following  equations: 


0.6  {2 


Rg  =  9.63x10  -A/^,  and  -  1.78x10  "A/, 


a.)  291 


f-O; 


l'J-24 


The  pre-exponent  in  the  expression  is  consistent  with  values  in  the  literature*’ (9  106  .s 
10'^  to  1.57  X  10'~)  for  PS^oluene  solutions,  although  there  seem  to  be  two  sets  of  values  clu:>- 
tered  around  1.0  x  10'~  and  1.5  x  10'“.  The  same  references  repon  a  range  of  exponent  values 
from  0.623  to  0.579,  which  is  consistent  with  the  result  reported  here  and  the  expected  theoreti¬ 
cal  value  of  0.6  (section  1,2. 1.3),  but  again  the  values  seem  to  clu.ster  around  two  values  of  0  58 
and  0.61. 


3.00- 


!og  Rg 


1.00 


-1.00 


log  A2  ,3.00 


-S.00 


5.00 


5.50 


5.00 

log  M, 


650 


7.00 


w 


logjQ  (circles)  and  log^Q  (squares)  as  a  function  of  1o2|q  .M  .  for  samples 
N170K,‘N770K,  N2M,  and  N3M. 


Figure  6.  Plot  of  log  R  and  log  A«  versus  log  M 

9  w 

The  A2  relationship  reported  here  is  also  in  agreement  w-ith  others"^  '  in  which  the  pre¬ 
exponential  factor  ranges  from  0.0158  to  0.0281.  and  the  corresponding  exponents  from  -0.286 
to  -0.329.  There  were  two  reponed  relationships*^'"'  with  significantly  smaller  pre-exponential 
factor  and  exponents;  A2  =  0.00636  and  A-,  =  0.00436  Both  repon  measure¬ 

ments  made  at  20°C,  which  would  correspond  to  smaller  A-,  values,"^  but  one  of  the  other  ref¬ 
erences^**  which  was  in  agreement  with  the  present  work  was  also  at  20°C,  so  there  doesn't  ap¬ 
pear  to  be  a  clear  correlation  with  temperature.  The  latter  two  expressions  come  closer  to  the 
theoretical  exponent  of  0.2,  predicted  by  the  two-parameter  theory,"^ 


3.2  Dynamic  Light  Scattering 

Dynamic  light  scattering  measurements  were  made  on  all  samples  except  N170K.  (For  the  la¬ 
ser  employed  in  this  study,  the  scattered  light  intensity  from  a  low  MW  sample  such  as  N170K 
is  too  low  to  obtain  a  good  correlation  function.)  A  typical  correlation  function  for  sample 
BR3M  is  shown  in  Figure  4.  The  DLS  data  were  obtained  for  rwo  purposes:  ( 1 )  to  extract 
infinite  dilution  values  of  the  diffusion  coefficients  (D^^)  for  unimodal  samples  so  that  kT  and  b 
values  from  Equation  9  could  be  evaluated,  which  are  needed  for  molecular  weight  distnbution 
transformation,  and  (2)  to  see  if  the  panicle  size  distrubution  (PSD)  data  analysis  routines  from 
Brookhaven  Instruments  (BIC)  could  differentiate  between  unimodal  narrow,  unimodal  broad, 
and  bimodal  distmbutions. 
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3.2. 1  Diffusion  Coefficients 


Autocorrelation  functions  measured  at  a  given  scanering  angle  and  solution  concentration  yield 
an  average  linewidth,  F.  which  can  be  converted  to  a  diffusion  coefficient  via  Equation  1 7  For 
monodisperse  systems,  this^ extrapolated  diffusion  coefficient  should  be  independent  of  the  scat¬ 
tering  angle,  through  the  q“  correction.  This  angular  independence  was  verified  experimentally 
for  sample  N3M,  and  thus  diffusion  coefficients  were  measured  at  one  or  two  angles,  luid  the 
values  averaged.  The  diffusion  coefficients  were  also  measured  as  a  function  of  concentration 
for  samples  BR2M  and  N3M,  and  values  of  kj^  (Equation  8)  were  determined  to  be  227  and  2^4 
mL/g,  respectively.  The  kj^  values  determined  in  this  study  are  consistent  with  literature  values 
for  PS/Toluene  solutions,  as  .shown  m  Figure  7.  There  doesn't  seem  to  be  any  theoretical  basis 

3.0 
2.5 
20 

log  kg  1.5 
1.0 
o.s 
0.0 

4.0  4.5  5.0  5.5  E.O  6.5  7.0 

log  M 

Plot  of  logjQ  kg  vs.  logjQ  for  this  work^  *  Ref  20  ||[| );  Ref  21  A  i: 

Ref  23  );  and  Ref.  27  {Ifr  ).  Dotted  line  is  ti,c  best  fit  to  all  data. 

Figure  7,  Plot  of  log  versus  log  M 

VJ  Vr 

for  plotting  the  log  of  molecular  weight  versus  the  log  of  kj^;  however,  a  reasonable  correlation 
is  obtained  which  allows  extrapolation  of  kjj  values  for  samples  in  which  a  concentration  de¬ 
pendence  was  not  measured. 

Values  of  the  infinite  dilution  diffusion  coefficient  (D^),  the  kp  value  used  to  conven  D(c)  to 
Dg,  the  hydronamic  radius  (Rj^)  determined  from  via  Equation  4,  and  p  =  Rg,^!,  are  pre¬ 
sented  in  Table  3. 

The  diffusion  coefficients  for  the  unimodal  (narrow  and  broad)  M’WD  polymers  were  deter¬ 
mined  from  a  single  exponential  fit  to  the  autocorrelation  function,  which  yielded  an  average  F, 
which  was  then  convened  to  D  via  D=Fq",  The  diffusion  coefficients  were  convened  to  infi¬ 
nite  dilution  values  using  either  a  measured  kg  value  or  a  value  extrapolated  from  Figure  7, 
based  on  the  determined  from  static  light  scattering. 

Tfre  analysis  of  the  bimodal  samples  was  more  indirect,  because  they  were  analyzed  using  the 
double  exponential  fitting  procedure,  which  yields  panicle  diameters,  not  linewidths.  Two  di¬ 
ameters  were  obtained  from  the  data  analysis,  and  were  convened  to  diffusion  coefficients  via 
the  Stokes-Emstein  equation.  The  two  diffusion  coefficients  for  each  bimodal  sample  were 
converted  to  values  using  a  single  kg  value  (from  Figure  7,  based  on  the  M^,  detemiined 
from  static  light  scattering)  and  the  total  solution  concentration.  Usinc  averase  M  and  kn, 
values  for  a  bimodal  system  is  a  rough  approximation,  however  an  alternate  approach  is  to  use 
kg  values  appropriate  for  each  component,  coupled  with  the  concentration  of  that  component. 
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This  analysis  yielded  similar  results.  Probably  the  most  realistic  approach  wouiii  he  to  caJcuiatc 
an  effective  concentration  of  each  component  (which  is  less  than  the  overall  concentration,  hut 
more  than  the  indivudual  concentration),  to  be  used  with  the  kp  value  associated  with  the  par¬ 
ticular  molecular  weisht  involved.  Needless  to  sav.  the  D  ^  and  the  Revalues  obtained  for  the 
^  'Oh 

bimodal  samples  are  very  approximate,  but  are  included  for  comparison. 


Table  3.  Dynamic  light  scattering  parameters 


■ . . . —  .  . .  . . . 

10^x  Do’ 

cm^sec"^ 

If 

Kq, 

3  -1 
cm  g 

nm 

(=  Rg'R,) 

N770K 

1.663 

105* 

23.9 

;  1.531 

BR2M 

0.9155 

227  (251*) 

43.35 

I  1,631 

N2M 

0.8698 

282* 

45.63 

I  1.637 

BR3M 

0.7104 

294** (372*) 

55.87 

1  1.640 

N3M 

0.7094 

294  (378*) 

55.94 

1.605 

2.874 

105* 

13.6 

1.173 

_ 1 

51.50 

1 ,373 

—  1 

1  1,987 

115.6 

1  1,582 

*estimated  via  Figure  7 
♦^measured  for  sample  N3M 


The  values  from  the  narrow  and  the  broad  distributions  are  plotted  against  molecular  weight 
in  a  log-log  plot  in  Figure  8.  The  two  broad  distribution  samples  appear  to  follow  the  linear 
relationship  of  the  narrow  distribution  polymers.  TTie  Rj^  values  from  Table  ?  can  be  plotted 
similarly  with  molecular  weight  (not  shown  here),  although  since  the  radii  are  calculated  di¬ 
rectly  from  the  diffusion  coefficients,  the  data  is  not  really  independent.  The  dynamic  relation¬ 
ships  derived  from  these  plots  are 

Do  =  3. 18  a:  .  and  =  7. 68  a:  10" (27) 


The  pre-exponential  factor  for  the  diffusion  relationship  reponed  here  is  a  bit  smaller  than  most 
of  the  literature  values  **  “*  which  range  from  3.4  -  3.7  x  10  ,  although  one  value""  is 
even  smaller  at  2.29  x  10“*.  The  absolute  value  of  the  exponent  in  this  diffusion  relationship  is 
also  smaller  than  those  reported  for  the  same  references  (0.577  -  0.587),  with  the  same  excep¬ 
tion  that  Seery  et  al.~^  reported  an  even  lower  value  of  0.533. 
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Two  papers  report  relationships  with  larger  pre-exponent lal  taciors  < 0  0 1 3  1  .md  0  ()  1  (*3  < 
and  smaller  exponents  (0.560  and  0.577,  respectively).  Tlie  two  effects  counteract  eacli  <>ther 
somewhat,  so  overall,  the  data  may  not  be  that  different. 


Figure  8.  Plot  of  log  versus  log 

The  p  data  in  Table  3  are  reasonably  constant  (average  p  =  1.603)  except  for  the  bimodal  sam¬ 
ples,  in  which  both  and  Rj^  are  very  approximate.  (The  R„  values  for  the  bimodal  samples 
were  determined  frortf  Equation  26,  using  either  a  deterrnined  from  static  light  scanering, 

or  the  manufacturer’s  M^.)  This  value  of  1.603  is  in  good  agreement  with  the  literature  ;  Park 
et  al?^  report  a  value  of^l.59.  while  1.6  and  1.5  can  be  deduced  from  the  data  of  Appelt  and 
Meyerhoff^^  and  Varma,"^  respectively.  The  theory  regarding  whether  p  should  be  a  constant 
value  for  all  polymers  in  good  solvents  is  unclear."'^ 

3.2.2  Particle  Size  Distribution  Analysis 

Autocorrelation  functions  measured  for  the  narrow  distribution  samples  N770K,  N2M  and 
N3M  were  analyzed  by  the  various  particle  size  distribution  (PSD)  analysis  routines  described 
in  section  1.3.1.  The  results  were  consistent  in  that  unimodal  narrow  distributions  were  ob¬ 
tained.  The  double-exponential  analysis  was  either  not  available  as  an  option  (the  experimental 
data  could  not  be  fit  to  a  double  exponential)  or  a  double  exponential  fit  would  yield  a  major 
peak  (-99%)  at  the  expected  particle  diameter  and  a  very  minor  peak  at  an  unrealistic  (either 
too  small  or  too  large)  diameter.  The  peak  diameter  and  the  breadth  of  the  distribution  were 
invariant  with  the  tj^  of  baseline,  the  number  of  data  points,  and  the  method  of  data  anlysis 
used.  A  typical  distribution  is  shown  in  Figure  9a. 

Similar  results  were  obtained  for  the  broad  distribution  samples  (see  Figure  9b).  although  there 
was  considerable  variation  in  the  breadth  of  the  distributions  determined.  The  data  analysis 
routines  do  not  appear  to  be  sensitive  enough  to  distinguish  between  a  narrow  and  a  moderately 
broad  distribution.  The  distributions  computed  for  a  narrow  and  a  broad  molecular  weight  stun- 
ple  of  approximately  the  same  molecular  weight  (shown  in  Figures  9a  and  9b,  respectively),  ap¬ 
pear  to  have  equal  breadth.  Due  to  the  low  power  laser  employed  in  this  study,  and  the  neces¬ 
sity  of  keeping  the  solution  concentrations  below  the  overlap  concentration,  c*.  the  scattering 
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intensity  for  all  correlation  functions  vv;ls  vety  low.  Data  had  to  be  accumulated  for  lone  peri¬ 
ods  of  time  (30-60  minutes)  to  obtam  smooth  correlation  functions.  Possibly,  smoother  data 
w'ould  produce  more  consistent  results,  but  m  general,  the  technique  of  DLS  is  not  suitable  for 
obtaining  detailed  information  about  size  distributions  (see.  for  example.  Ref.  5.  10-14), 


G(d) 


G(d) 


logjQ  diameter,  nm 


G(d) 


Gtd) 


logjQ  diameter,  nm 


Particle  size  distributions  G(d)  versus  JogjQ  particle  diameters  (in  nm)  for;  (a) 
N3M  analyzed  by  the  non-negatively  constrained  least  squares  (NNLS)  routine; 
(b)  BR3M  analyzed  by  NNLS;  (c)  BI850K  analyzed  by  the  exponential  sam¬ 
pling  (EXPSAM)  technique;  and  (d)  BI4M  analyzed  by  EXPSAM. 


Figure  9.  Particle  size  distribution  results 


The  analysis  of  the  bimodal  samples  always  yielded  a  bimodal  distribution  w  ith  the  peaks  at 
their  approximately  expected  locations.  These  correlation  functions  were  analyzed  by  different 
data  routines,  and  essentially  the  same  results  were  obtained.  .All  distributions  shown  in  Figure 
9  were  analyzed  via  the  exponential  sampling  technique  or  the  non-negatively  constrained  least 
squares  analysis,  because  the  other  data  analysis  techniques  mentioned  in  section  1.3.1  (cumu- 
lants  and  double  exponential)  yield  only  one  or  two  size  values,  not  a  distribution,  however  the 
numbers  obtained  with  the  other  data  analyses  were  consistent  with  those  shown  in  the  figure. 

The  different  data  analysis  routines  appear  to  produce  distributions  with  their  own  characteristic 
shapes  (notice  the  similarities  between  Figures  9a  and  9b,  and  between  Figures  9c  and  9d).  The 
Exponential  Sampling  (EXPSAM)  method  tended  to  produce  smoother  bell-shaped  distribu¬ 
tions  than  the  non-negatively  constrained  least  squares  (NNL)  analyses,  which  had  a  more 
pointed  shape.  In  the  case  of  the  bimodal  distributions,  the  NNL  routines  always  produced  two 
distinct  peaks  which  did  not  overlap,  while  the  EXPSAM  results  were  more  continuous  over  the 
entire  size  distribution,  as  shown  in  Figures  9c  and  9d. 
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3.2.3  Molecular  Weight  Distributions 

Molecular  weight  distributions  (MWD's)  were  determined  for  the  four  polymer  s.miples  de¬ 
picted  in  Figure  9  by  first  transforming  the  PSD’s  to  distributions  of  linewitiths  ‘.r\}  via  the 
Stokes-Einstein  equation  (eq.  4)  and  then  transformuig  the  linewidth  distributions  to  MWD’s 
via  the  equations  given  in  section  1.2. 1.3.  AU  the  MWD  transformations  were  perfonned  on 
PSD’s  determined  via  the  NNL  analysis  routme,  because  the  MWD's  produced  from  the  EXP- 
SAM  technique  were  not  smooth. 

f 

Figure  10  show's  the  results  of  the  molecular  weight  transformations.  Tlie  values  of  the  parame¬ 
ters  kp.  ky.  ky’.  b  and  b'  affect  the  shape  of  the  resultant  distribution  as  well  as  the  numbers 
derived  from  it  (radius  of  gyration  and  molecular  weight  averages).  As  i  first  attempt  at  a  ® 

transformation,  values  of  kj^  were  estimated  and  the  experimental  values  of  the  other  four  pa¬ 
rameters  (as  reported  in  Equations  26  and  27)  were  used.  These  a  priori  values  prrxluced  a  re¬ 
alistic  MWD  for  only  one  of  the  particle  size  dismibutions  (N3M).  All  the  other  distributions 
required  some  parameter  adjustment  to  produce  a  smooth  distribution  with  meaningful  num¬ 
bers.  A  computer  program  was  written  to  calcualte  a  MWD  and  the  resultant  averages  fR^,. 

M^,  Mj^)  for  a  given  set  of  adjustable  parameters.  ~ 


log, Q  MW  log,,,  MW 


Molecular  weight  distributin  f(M)  versus  log,y  molecular  weight  for  .  (a) 

N3M  ;  (b)  BR3M;  (c)  BI850K;  and  <d)  BI4M.  All  were  rrans formed  from 
PSD’s  computed  via  the  NNL  data  analysis  routine. 

Figure  10.  Molecular  weight  distribution  results 

The  dynamic  parameters  ky  and  b,  transform  the  diffusion  coefficients  into  molecular  weight 
values,  and  thus  have  a  significant  effect  on  the  distribution.  For  the  bimodal  siunples.  increas¬ 
ing  b  from  the  measured  value  of  0.563  to  0.58  (which  is  more  consistent  with  the  literature) 
produced  a  smoother  distriubtion.  The  B1850K  sample  required  a  furthc  minor  adjustment  of 
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parameter  ky  to  3.1  x  10"^  (from  3,18  x  Ur'^i  to  produce  the  ni.al  ikstnl’iuiion  Fhe  ^!a^c  pvi- 
rruTieters  k.j-'  and  b’,  which  transfomt  molecular  weiehi  valucN  into  radn  of  evranori.  *1111%  af¬ 
fected  the  final  values,  so  they  were  left  at  their  expenmentailv  deiennined  values 

C 

The  second  virial  coefficient  for  diffjsion.  kj^.  transfomrs  the  diffusion  coefficients  measured 
at  finite  concentrations  to  inifnite  dilution  values.  D^.  Of  the  four  distributions  shown  in  Fips 
9  and  10.  kjj  w'as  only  measured  for  sample  N3M.  and  this  value  of  kj^  w  as  used  unchanged  m 
the  molecular  weight  transformation  process.  The  other  k^  values  were  detemnned  by  trial  and 
error  using  a  value  estimated  from  Figure  7  as  a  startmg  point.  Tlie  final  values  were  .ill  less 
than  those  estimated  from  the  figure.  For  sample  BR3M.  the  value  of  kj^  was  decreased  h\ 
about  30%.  and  for  the  two  bimodal  polymers,  approximately  15%  each  !Mcident.iily ,  tiic 
measured  value  of  kj^  for  sam.ple  N3M  is  about  30%  lower  than  the  value  calculated  from 
Figure  7. 

Table  4  contains  the  molecular  weights  and  radii  of  gyration  calculated  from  tire  distnbutions 
shown  in  Figure  10,  The  numbers  calculated  from  the  unimodal  distributions  matched  the  ex¬ 
perimental  data  better  than  the  two  bimodal  samples.  .An  attempt  was  made  to  use  tw  o  different 
kp  values  in  the  molecular  weight  transformation  process  for  the  bimodal  samples,  but  there 
W'as  little  effect  on  the  final  distribution. 

Table  4.  Parameters  calculated  from  molecular  weight  distributions 


HSRH 

peak  MW’s 
(X  10'®) 

Adjusted 
'  Parametes 

N3M 

2.99  (+3) 

1.14 

86.7  (-3) 

‘  none 

BR3M 

2.84  (-10) 

1.11 

88.1  (-4) 

2.25 

kp  =>  250 

BI850K 

1.21  (+44) 

2.96 

5^.0  (-13) 

0.161  (-4). 
1.620  (-24) 

:  kp  =5  b() 
b  =5  0.58 
kj  =0  3.1E-4 

BI4M 

5.16  (+26) 

3.43 

128.3  (-  1) 

0.437  (-56). 
6.37  (-32) 

;  kp  =>  350 
b  =5  0.58 

Numbers  in  parentheses  indicate  percent  differences  from  expiected  values 


One  interesting  thing  to  note  for  both  bimodal  samples  is  that  even  though  the  calculated  over¬ 
all  weight -average  molecular  weights  are  significantly  higher  than  expected,  the  peak  molecuku 
weights  for  each  of  the  components  are  lower  than  exp>ected.  This  unplies  that  either  the  tw  o 
molecular  weight  fractions  are  broad  di.  tribution  with  a  high  molecular  weight  tail,  or  that  the 
estimate  of  the  fraction  of  the  larger  molecular  weight  component  is  too  high.  Estimates  of  the 
polydispersity  of  each  of  the  fractions  are  about  1.05,  thus  ruling  out  the  first  possibil¬ 

ity.  The  mole  fraction  of  each  component  can  be  determined  by  dividing  the  sum  of  the  fi  M^>'s 
for  each  component  by  the  sum  of  the  f(M.)‘s  for  the  entire  distribution.  Using  this  approach, 
the  mole  fraction  of  the  lower  molecular  weight  component  is  81%  for  the  Br4M  sample,  luid 
82%  for  the  BI850K  sample.  Converting  the  original  w'eight  fractions  of  each  component  (see 
Table  1)  to  mole  fractions,  the  lower  molecular  weight  components  should  be  94%  .md  bh'  i . 
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respectively  Tlie  cau'e  ot  this  skcutiie  ctt  the  ilistnhuia'ii  tovs.ud  tlic  luchei  iiuiicvihaj  ■Anyiit 
components  is  probably  caused  by  the  inherent  buis  ot  hcht  scattennc  unsajii  the  iuchc!  rnv 
leculor  weight  species 

4.  Conclusions 

The  static  light  scattering  experunents  were  straighitorward.  .md  .ue  leasouahh  vi»nsiv!cns  with 
literature  vales.  Difterences  in  the  refractive  mde.v  merernen!.  .uui  ihc  Rasleigii  ratio  ot  tlic 
calibration  liquid  used  in  the  data  analysis  could  account  foi  the  diftercnccs  nicse  saiue-  aic 
often  not  published  w  ith  the  light  scattering  data. 

Due  to  the  large  variability  in  the  published  DLS  data  for  toluene  solutions  ot  poKsiMenc.  it  ss 
clear  that  the  technique  of  DLS  for  even  a  sunple  analysis  of  diffusion  ccKtficients  is  not  a  nia 
hire  field,  and  expenmental  data  do  not  always  correlate  with  theoretical  predictions  Hie  rea 
son  for  the  large  discrepancy  m  published  data  ts  probably  due  to  the  fact  that  the  d>  nanne  light 
scattering  e.xperiment  u  not  easy  from  an  experimental  point  of  view  ,  or  from  the  daui  intcrpre 
tation  side. 

For  example,  the  selection  of  the  delay  time  <At)  to  be  used  when  collecting  data  tvn  the 
autocorrelation  function  affects  the  shape  of  the  correlation  function,  and  thus  the  munbeis  that 
are  derived  from  the  data  analysis.  For  systems  with  large  polydispersity,  the  use  of  muitipie 
sample  times  is  recommended  to  incorporate  all  possible  relaxations,  panicularly  when  the 
number  of  hardware  channels  (N)  is  limited.  (Multiple  sample  times  were  used  for  the  analysts 
of  sample  BI4M.)  If  multiple  sample  times  are  used,  and  too  large  a  range  of  sample  tunes 
incorporated,  information  on  the  larger  particles  can  be  lost,  whereas  if  the  range  of  sample 
tunes  is  too  short,  the  information  on  the  smaller  particles  is  sacrificed  Another  experimental 
parameter  that  affects  the  results  of  a  DLS  experunent  is  the  number  of  samples  used  to  obtain 
the  correlation  function.  /  c  ,  the  length  of  time  that  data  is  collected.  For  systern.s  w  ith  high 
scattering  power,  the  experiment  does  not  have  to  be  run  as  long  to  accumulate  the  s.mic  num¬ 
ber  of  samples  as  from  an  expennient  in  which  the  scattering  is  not  as  great  Ohviousls  .  the 
more  samples  that  are  accumulated  in  a  given  experiment,  the  smoother  the  correlation  function 
will  be.  and  the  more  preci.se  the  data  generated  from  that  experiment  will  be.  however  a  practi¬ 
cal  time  limit  for  the  data  collection  must  be  set.  Sometimes  a  shoner  experiment  duration  is 
desirable;  for  instance,  the  longer  the  experiment  is  run.  the  greater  the  chance  that  a  spurious 
dust  particle  will  diffuse  into  the  scattering  center,  or  that  the  ernironmental  coiulinons  (tem¬ 
perature,  humidity.  la.ser  intensity,  etc.)  wall  change. 

Section  1.3,1  bnefly  describes  some  data  analysis  routines  which  are  commonly  in  use  for  ana¬ 
lyzing  correlation  functions.  Many  other  methods  have  been  used,  and  are  currently  under  de 
velopment  -  in  fact,  much  of  the  current  DLS  literature  is  devoted  to  developing  new  mathe¬ 
matical  treatments  of  this  data,  and  of  assessing  the  accuracy  of  current  methods  Thus,  the 
analysts  of  DLS  data  is  still  an  area  of  active  research. 


In  conclusion.  DLS  is  well  suited  to  measure  the  dynamic  pro[x?nies  ot  a  moiiodispersc  pois 
mer  solution,  such  as  translational  diffusion,  and  hydrodynamic  radii  which  w  hen  coupled  w  uh 
information  from  static  light  scattermg,  provides  a  complete  dilute  solution  analysis  o!  a  poK- 
mer  in  a  given  solvent.  Indeed,  several  researchers"*^  "”  are  developing  combined  static  and  d% 
namic  light  scattermg  spectrometers  so  that  the  static  and  dynamic  paiaineiers  may  !>e  obtained 
from  the  same  expieriment. 

Conversion  of  an  experimentally  detemiined  Imewidth  distribution  to  a  distribution  ol  moiccu- 
lar  weights  is  possible  if  the  parameters  kj^.  Lp.  Lp’.  b  and  b'  are  known  or  can  he  estimated, 
and  accurate  values  of  M  ,  and  R  ,  are  available  from  sfaiic  Intht  scatterins  measurements  I’s- 
ing  DLS  to  extract  mformation  on  the  polydispcrsity  of  a  polymer  -  ample  should  be  limned  to 
determination  of  a  mean  value,  and  possibly  one  or  two  moments  about  that  mean,  due  to  the 
inherent  difficulty  of  extracting  a  unique  Imewidth  distnbution  from  an  experimental  correla¬ 
tion  function.  Unimodal  distnbutions  can  be  distmguished  from  bimodal  distributions  if  the 
two  species  are  at  least  a  decade  apart  (more  closely  spaced  distributions  may  be  distinguish¬ 
able  under  certain  conditions),  but  extraction  of  detailed  mfomiaiion  aboui  the  individual  frac¬ 
tions  is  probably  unreliable. 
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